


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1969 


The cross section for the formation of H 
[Superscript +] [Subscript 2] in the reaction of 
fast protons with methane 


Smyth, Norman Robert Anton 


Monterey, California. U.S. Naval Postgraduate School 


http://hdl.handle.net/10945/11903 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| г D U DLE Y research materials and institutional publications created by the NPS community. 
FW узу, Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


ІШ KNOX appointed — and published — scholarly author. 

OM LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


NPS ARCHIVE 
1969 
SMYTH, N. 


THE CROSS SECTION FOR THE FORMATION OF Ht 
IN THE REACTION OF FAST PROTONS 
WITH METHANE 
by 


Norman Robert Anton Smyth 











DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CA 93943-5104 


United States 
Naval Postgraduate School 





THESIS 


THE CROSS SECTION FOR THE FORMATION OF Н, 


IN THE REACTION OF FAST PROTONS WITH METHANE 


by 


Norman Robert Anton Smyth 





June 1969 


This document has been approved for public re- 
Lease and sale; its distribution is unlimited., 


Le 1 
NAVM, PUSTZRAATULUS и. 
Бакру САР, 27840 





8 + 
The Cross Section for the"pemmetion of Нн, 


іп the Reaction of Fast Protons with Methane 


by 


Norman Robert Anton Smyth 
Captain, Canadian Armed Forces 
B.Sc., University of Alberta, 1963 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN PHYSICS 


from the 


NAVAL POSTGRADUATE SCHOOL 
June 1969 


ABSTRACT 


The capture cross section for the formation of H, in the 
reaction H' + CH, 7 H, + CH. was measured at incident proton 
energies of 70, 85, 100, 150 and 200 eV and covering the scat- 
tering angles of 5° to 49.59 (lab coordinates). At 100 eV and 
below the curve of the cross section versus angle shows a sharp 
peak at about 46° whose position approaches the theoretical limit 
of 46.9° with increasing energy. Above 100 eV the peak was too 
small to be observed and only an upper limit can be placed on 
the value of the cross section. Typical values of the total 
cross section are 2.0 x 105 ~ at 70. ev anadai Ork 19758 at 
100 eV. The magnitude and energy dependence of the cross section 
as well as the angular position of the peak all are in essential 


agreement with the classical theory of ion-molecule rearrangement 


collisions proposed by Bates, Cook and Smith. 
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I. INTRODUCTION 


The correct description of complex atomic and molecular inter- 
actions is one of the major achievements of modern quantum theory. 
While there appears to be little doubt about the correctness of 
the theory in principle, it has become evident in practice that 
the mathematical complexity is such that only the very simplest 
cases can be treated with anything approaching complete rigor. 
Therefore many approximations are employed to simplify these 
mathematical computations. The question then, of the accuracy, 
nature and range of applicability of these approximations is of 
considerable practical importance. In general these questions can 
only be answered by comparison of the experimental results and the 
predictionsof a specific approximate calculation. 

This paper presents the results of a series of measurements 
which confirm the predictions of a classical theory of ion- 
molecule rearrangement collisions at high impact energies as pro- 
posed by Bates, Cook and Smith [1]. This theory uses an impulse 
type approximation to obtain the cross section for the capture of 
a light atom or ion from a target molecule by a fast moving pro- 
jectile. Applying this theory to the reaction: 


H + CH, Hj * CH (1) 
н + Сн, = Ho + CH. 


the cross section is predicted to show a sharp peak at 46.9° with 


an upper limit on the magnitude of the cross section of 


-20 
С 


2 4 е à 
П... Грос uh m at 100 eV incident proton energy. It is further 


11 


predicted that the magnitude of the cross section should decrease 
rapidly with increasing proton energy, approaching an energy de- 
pendence of в 2-2 asymptomatically at high energies (above 500 еу). 
The range of validity of the theory extends from about’ 50° to oU0me® 
when applied to reaction (1). 

Our measurements were carried out at 70, 85, 100, 150, and 
200 eV, covering the angles of scatter from ж to Е 7 We 
found essential agreement between the predictions of the theory 
and the experimental data, indicating that the classical impulse 


approximation proposed by Bates, Cook and Smith is a valid model 


for reaction (1) in the energy region considered. 


11% 


LL THEORY 


A. CLASSICAL THEORY OF ION-MOLECULE REARRANGEMENT COLLISIONS 

A classical theory of ion-molecule rearrangement collisions 
at high impact energies has been proposed by Bates, Cook and 
Smith [1]. The impact energies are assumed to be high enough 
so polarization forces and chemical binding energies of the 
colliding molecules can be ignored. The theory applies to various 
ion molecule rearrangement reactions of the type: 


t fo ! T t 
X'X. * YZ Z 8012277 ШІ! 





мһеге Х. апа Ү аге simple atoms or ions. The bar indicates which 
ions are fast in the laboratory coordinate system. An example of 
a rearrangement reaction of this type is: 


Hi + сн, 5 ЧО ссн (2) 
EOD Kee 


The basic assumption made by Bates, Cook and Smith, is that 
process (1) may be described by a classical impulse approximation 
similar to that developed by Thomas [2], for the description of 
electron capture, where each composite system is regarded as a 
И + + 

loose cluster of atoms and ions. (Eg H, = H+ H anad 
CH, = С + ЗН + Н). 

According to this model, we can consider the reaction to occur 


in the following manner: consider a particle of mass M, moving with 


1 


velocity V1 through molecules of loosely bound particles of mass 


M, and Мк. Figure 1 shows the sequence of collisions that lead to 


TS 





FIGURE | THE CAPTURE MECHANISM 
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First mass M. collides with mass M.. 


the capture of М„ by Mi. 1 2 


2 


1 


This causes М. to be scattered at an angle 91 and velocity Va 


while M, recoils at an angle 9, with velocity VA. Mass M, then 


suffers a second binary collision with mass M This causes M 


к” 2 
to scatter at an angle 95 with velocity LE 
Now if: (а) Of = Of + 01 (3) 
2 2 1 E. Б 
Or vi = vy! 


ке 22 oN 
And if: (b) уу = V5 


then masses М, апа M, may have a relative energy of motion below 


that required for separation, and so, the two particles may com- 
bine and proceed on as one. If the particle with mass М, 15 а 


composite particle then its disruption can be avoided if 


> = 

Е (4) 
vi ~ КТ . 
In the particular case of: 


М2 М) << М, 


condition (3) is satisfied if 9. = 4,59 and 95 м 90? (see Section 
II of this Chapter), but then condition (4) is violated. In the 


case of the reaction: 


н? ТЕСТКЕ H =) CH 
D 4 —2 3 
rY = = = = * ° в 
мһеге М. М M, M, and M Мен ‚ һе violation of condition 


2 


(4) is not relevant since M, is not a composite particle. 


1 


Consider now, the first binary collision between M, and M, às 


shown in Fig. 2. The probability of scattering M, into the solid 


-— > -— 
angle ай(04) with velocity between уд апа V. t dvs is: 


а = 9,4(92) d0(07) (5) 


15 





FIGURE 2. THE FIRST BINARY COLLISION 





FIGURE 3. THE SECOND BINARY COLLISION 
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where 012192) is the appropriate differential scattering cross 

section. Since the solid angle is the cone of semi-angle 9, within 
! . 

de; then: 


jo res : t t 
dQ(9;) 27 sin ө; do 5 


and equation (5) becomes: 


- ! : ! t 
q 270,,(95) sin 0, dos . (6) 


Consider now, the second binary collision between M, and M, as 


shown in Fig. 3. The probability of M, having an impact parameter 


2 
between @ and@+ d@ and an azimuthal angle between Фф апа №) + аф 


at a distance r from M, is given by: 


И ар 
2 5 249 - --- (7) 
mer 
where r is the М-М, equilibrium internuclear separation. But 
particles incident on М, with impact parameter @ and azimuthal 
angle №} аке scattered into the solid angle dr(05) at O75. 
Hence: бае = с (95) 40(65) 
where АС is the appropriate differential scattering cross 
section. 
° @ "у = 1 11 " 
But: а (95) sin 9 de, 


Hence equation (7) becomes: 


| | 
раж о кси ---; (8) 
Дпх 
Now, for М, to capture M, ; Vi - v5 the relative velocity of Mi 


and M, must be contained within a volume of velocity space 


Я 


determined by D, the mutual affinity of these two systems. Іп 
the high velocity 1ітіту (Ва се, СООК апа БТЕ ЕТИ 


requirement as: 


Da 
2 в. 2D 
(v5) sin 05 doy ch) ы = EU v) (9) 
where ха = М. in the case of МЕ (10) 


Substituting equation (9) into equation (8) we get 


се (Ө) 372 
DEO А (11) 


2 
" " 
3r (v5) dv5 


Now, the capture cross section arising from the M, -M,-M, 


sequence of binary collisions is: 


СЕМ ВМ Е (12) 


where У 15 a dimensionless factor, less than, or equal to unity, 
which allows for the fact that M4 and М, пау approach each other 
such that in the resulting state, the affinity is less than that 


assumed. 


Substituting equations (6) and (11) into (12) one gets: 


= = = ИИИ ЕЕ -v 1 T po e mm 
Q(M,-M,-M,) = YA 5 віп 0) 40, "mo дот. (13) 


ло a) op M? 95 
ar 


Changing from Lab. to Center of Mass coordinates, one gets for 


the case of M., = М. << М 





1 2 k 
ТС” 3/2 
O(M.-M.-M.) = ¥ 12 2k 2D 14 
ЕМ = 3r? - D 
2 1 
where c and c denote the differential scattering cross 
sections in the center of mass coordinate systems. This 


18 


differential cross section is at an energy of relative motion of 
м м, and the affinity D is at an М.-М, internuclear distance 

оф М2 т. To obtain the total cross section for capture, Q, one 
now sums over all binary collision sequences which lead to the 
same final capture process. For example, the total cross section 


for process (2) is 


Q = 4Q(H-H-c) + LAH H C). 


This represents the two possible collision sequences shown below: 





t H 
H + 
H 
H 
He 
H HO 
Sc 
HO но 
OH OH 
+ + 
The (H - H - C) Sequence Тһе (Н- Н - C) Sequence 


Solution of Rutherford's expression for the differential 


scattering cross section when v, is high enough for the atomic 


1 
field to be Couloumbic, shows that: 


-11 
Q(M,-Mo-M.) vi . (15) 


The range of validity of (14) does not extend indefinitely 


I 
(i) The De. Broglie wavelength limitations, 


as v, increases. Two requirements must be satisfied: 


(ii) The Heisenberg uncertainty relations, 


19 


Bates, Cook and Smith [1], computed the required differential 
scattering cross sections о ,(907) and о „ (907) using classical 
methods and the analytic representation of the relevant Hartree 
potentials as given by Byatt [3]. From the results of these 
computations listed in Ref. 1, the cross section for the re- 
arrangement reaction: 

H ОСН н> ЕСН 
З с 


is determined to be as shown in Fig. 4. 


- - 
В. НЕР emer m H, T = REACTION KINETICS 


The rearrangement mechanism for the reaction is shown in 


d + 
Fig. 1, where Mi represents the H , М> the H and М represents 


the CH. To determine the exact value of the scattering angle Ө! 


3 1 


one has to consider the kinematics of the arrangement process. 
If the speed У is assumed to be large enough so the binding 


energy of M, to M, can be ignored, then the conservation laws 


k 


can be applied in the usual manner. 


Consider the first binary collision of the reaction shown in 


Fig. 1. From conservation of energy and momentum we have: 
муу = MICA 39S CD (16) 
Mv, = МУ] cos 91 ~ MV cos ө; (17) 
О = МУ sin Ө] - MV sin 95 (18) 
Now: Mi = М> = М 
Hence: the above equations become 
У = (01) + (02) (19) 


20 


(04,0 Х б стг 





47 


-8 


ҒКОМ 
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H* 


FIGURE 4 
PREDICTED CROSS SECTION FOR 


+ CH, 
BATES 


+ 
— Н, + CH, 


COOK SMITH 


3 Log Energy (м) 


О cos O "РЕЈЕСО5 0! (20) 


1 1 1 2 2 
= 1 1 yc 1 1 1 
O v} Sin Ө) v) Sin 9> (20) 


Now ө; can be eliminated between (20) and (21) by rearranging, 


Squaring, and adding the two equations to yield: 


2 1 2 1 т = 1 2 
У + (уу) 2У1У4 cos 9, = (v5) 
2 2 2 
e T = = 1 
But from (19): (v5) V1 (уу) 
2 2 2 2 
. + 1 28; 1 - іг 1 
Thus: Д (уу) 2v. VY cos e: V1 2) 
2 
: 1 = ! т 
Непсе: (v,) У1У1 Cos 9, 
ыды vi = у, соз 01 (7224) 


Substituting equation (22) back into equation (19) we get: 


а 2 A de 
un Д (cos 91) + (v5) 
2 2 E 
e — 1 = 1 
On. a (1-cos 91) (v5) 
Я t = 1 1 
V5 v, Sin 9, (23) 


Equations (22) ана (23) Give Us vi and у» in terms of Vy and 


өү. Now we wish to obtain V1 and V5 in terms of ul and 95. Thus 


we arrange equation (20) and (21) such that when we square and 


add them, we eliminate e: and so obtain: 


2 2 2 
t " 1 ' = 1 
QM B (v5) 20175 cos ө} оа 


po 


2 2 2 


But from (19) we have: (v1) = З= (v5) 
° 2 1 2 1 ' = 2 1 2 
Thus: es (v5) - 2V4 V5 cos 9, Vi ( v5) 
к I LR 
Hence: ш ас (24) 


Substituting equation (24) back into equation (19) gives: 


v E CM T v cos? 9, 
OT: v; (1-cos*03) = (1) 
2% СЕЕ cee (25) 
Consider now, the second binary collision. Again ignoring the 


binding energies we can write the conservation of energy and 


momentum equations as: 


2 2 2 
= Мау.) = Evo) + M, Vy (26) 
a " " 
М.У МУМ" СОЗ 9% + М Му СОЗ 9, (27) 
= " 1 ML 1 
D Nun CMM cha (28) 
Now: М, = М, z М, 
а; - 
an Mi, Meg 
3 
Thus the above three equations reduce to: 
M 
2 2 k 2 
1 = 11 - — 
ie ME EI v (29) 
H 
Мк 
1 == 11 11 ee 
( v5) V5 cos e * M. Ук cos Ө, (30) 


23 


= 


k : 
O v5 sin о; у Sin 9. i (og 


T 


As we are not interested in e. we Can eliminate it by 


suitably rearranging equations (30) and (31) and then squaring 


and adding. Thus we get: 


2 
M 
2 2 k 2 
! $ " _ 1 " On e ; 
(v3) (v5) 2v3 v," сов 95 M. Vk (32) 
But from equation (29): 
M 
k 2 2 
—= V = ( 9 КР (м i) 
k 2 
ү 
Thus (32) becomes: 
M M 
2 = k 2 k 2 
1 " E ту! и = ЕЕ 1 M " 
( v5) * (v5) 225 cos Ө " ( v5) Т (М) 
Н Н 
i 2 МЕ 2 Ma 
f D 4 " +“ = = "ы" о" = 0 
(v5) d M. (v5) 1 M урур соз Ө» 


2 
; . a 1 В 
Multiplying through by: М +ъх [us yields: 
k H 


КМ = 





2 
М..-М E Мы М 
HE +( =] - 2 = [4 | cos on = 0 (33) 
H k 29 2H k 
M M, -M 
Now, let: а = and b = 
+ + 
М Мы M, +М 
Then (33) becomes: 
uy : У» 
= - 21-т |а coso EERDER (34) 
V V5 2 


24 


But now, for capture, we require that: Ө" ы Ө! + Ө 








1 2 
. Has 1 1 
thus: cos 95 x COS (91 + 95) 
le., COS 05 - Cos 01 cos о, = Sin OF sin 95 and using (24) and 
(25) to eliminate 9; we have: 
t 1 
vo | al 
cos 0! ~ cos OF; — |- sin OF | — |. 
2 IN V I UV 
i ДІ 
Hence equation (34) becomes: 
V5 = yoo vil 
eT - 2а = = cos 0! - — -=> sin 0! sub = 0 
2 2 ү 2 M 
vuy! vity t 
= И 
or (v5) - 2ау> cos 91 = sin 91 ~ b(v1)^ - 0 ( 35) 


But now, the second capture criteria is that: У ~ vl 


But, from equation (22) we have: V1 = v} cos 91. 


Hence, we require for capture: V5 w V, Cos Of. ( 36) 


Substituting equations (22) and (36) into (35) we get: 


Q! = ! T к= 1 e Ө! = t = • 
Т COS 1 га; БҮР cos 9. лп 1 COS 11 b( v2) = 0 


Using (23) to eliminate У» from this equation we get: 


2 2 2 ; 2 2 2 
Or : t : t t О! | Т pu 
v,cos 91 2av,sinOj; [ v, sinOlcos oF v,Sin9, cos M bv,sin 94 О 


which reduces to: ve oe e: - b v "E 94 = 0s 


Hence, for capture to occur we must have: 


M +M 
| k H 
{ап Ө! == = — . (37) 
1 b M, - M, 


25 


Мои М = М = 15.03506 and M, = a800797eamu™" Thus thesamgle at 
k = Н 


which capture occurs is: 


Өт = 46.926" (38) 


Hence we expect the capture cross section to show a pronounced 
peak at this angle. From Fig. 4, the magnitude of the cross 
20 


section is predicted to be approximately 1.4 x 10. ona if the 


energy of the incident proton beam is 100 eV. 


C. THE RANGE OF VALIDITY OF THE CLASSICAL DESCRIPTION 

A lower limit on the range of validity of the classical 
treatment of the rearrangement collision is imposed by the 
assumption made in the theoretical development that the energies 
involved in the collision are much larger than the binding 
energies. Hence, when the energy of the incident proton beam Ei 
is less than, or equal to about 10 times the binding energy we 
expect the peak in the cross section to become less pronounced 
and more spread out. Thus the theory is valid when Ei = 50 eV. 

An upper limit on the range of validity of the theory is 
determined by the two requirements: 

(i) The Di Broghlie Wavelength A of any particle involved 
must be much less than the smallest distance involved 
in the collisioni. 

Hence we require: A << § (39) 
(ii) The uncertainty in energy introduced by specifying the 


transverse position of the proton must be much less than 


the binding energy of the final product. Now, from ( 39) 


26 


we know the uncertainty in the proton position A is less 
than the distance of closest approach S. 


But we can write the uncertainty relation: 


ћ 
ROADS 2 ( 40) 
as $Ар >> Ші 
2 
>> Й 
Thus Ap 25 (41) 
2 2 
ы жн E (42) 
H 8м,5 


We then demand that the binding energy D be greater than AE, or: 


(43) 





To find what restrictions equations (39) and (43) place on 
the energy of the incident proton beam we have to relate the 
distance of closest approach S to the energy of the proton beam 
E,. This can be done quite readily in the case of coulomb 


i} 


collisions. By considering the collisions to be coulomb it will 


be shown that: ee SN 
2m 
(> 
апа А = 
E: 


where Ci and с, are constants. Because S decreases faster than Л 


as Ei increases we see that equation (39) will indeed place an 


upper limit on the energy Е. 


27 


In a Coulomb collision the point of closest approach is 


yielded by solving the equation: 


2 
b 
1 - > -AS -o (44) 
U 
S Hi 
(see equation 3-4-3 of Ref. 4) 
— a 2 
04, the Center of Mass energy, can be expressed as Ui = 5 0 V eC 
where ш 15 the reduced mass and ons is the relative speed of 


approach. (S), the potential energy at the distance of closest 
approach is; 
_ (Ze ге! 
gis) = ZEZ e) 


S 


The impact parameter b is given by: 


1 1 (Н) 
LL 
"rel 
(see equation 3-8-8 of Ref. 4) 
Where (Н) 1s the Center of Mass scattering angle. 


Now we have two collisions to consider: 


25 : 
(i) In the first binary collision (H -H) we have (see Fig. 1) 


(Н) = 90°, пе амы = (46) 


~ " = ч — = 1 
(Н) з Өн = 900, ШЕМ уу M (47) 
and since Е. ~ ЬЕ, ме have у! ~ ES V 
p des] 2 (2 1 Е 
Making these substitutions in (45) we get: 
2 
E Б. 
H - H 1 
2 
b = 12 — (49) 
Н-С Е 
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where: Е = 22 M. V 


+ 
Hence for the H - H collision, equation (45) yields 


2 4 
саг тст но 

E, E^ 
i 

2 

Р І T e 
Which gives: S(H - H) ~ 2.4 rmm ( 50) 
1 


Similarily, for the H-C collision we get: 


<2 _ 1225 _ 144 ot A6 
E 2 
1 E, 
Е 
which gives: S(H-C) ~ 7.4 Е ° (51) 
1 


+ , . 

Hence the H - H collision has the smallest distance of closest 
approach and hence equation (50) will establish the upper limits 
of validity of the classical approximation. 


Consider first the restrictions imposed by equation (39). We 


require that: jl ES 
2 
E n << т 
2 
М, Ei l 
EG 
: Ei с 202 D Ergs (52) 
= 
Or E, CINES X 10 4 ave 


ў 


Thus the wavelength considerations require BE, << 1.45 x 10° ev if 


p 


our classical approximation is to be valid. 
Consider now the restrictions imposed by equation (43). We 


require that: „2 


р >> 
2 
8 

My 5 
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2 
А 
21 


Ог р >> ——. 
вм (2.4) e" 


ET 
The binding energy of H, is 2.6 eV. Thus we get: 


А 
2.4e 
Е, << EN D E ШЕСІ 


= 3.93ex me Ergs 


= 2 пе и wy. 


Hence equations ( 39) and (43) are satisfied if the energy of 
the incident proton beam is less than 1,000 electron volts. Thus 
the classical description of process (2) is valid in the 75 to 


200 eV range where the experiment was performed. 
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ТТТ. EXPERIMENTAL APPARATUS 


The experimental apparatus used in this experiment is the 
result of the work of many people. The overall apparatus, as 
shown in Fig. 5, was assembled by Bush [5], who demonstrated its 
Satisfactory operation. 

A reasonably mono-energetic beam of hydrogen ions (AE — 2 eV) 
is produced in the Duoplasmatron, which was built and investigated 
by Carter [6]. The beam is mass analyzed in the mass spectrometer 
constructed by Strohshal [7], and н” ions proceed on to the scat- 
tering chamber which contains the target gas, c The scattering 
cell is placed on the magnetic axis of a large cylindrical magnet 
at the position of maximum field. H, ions, produced by the col- 
Ето ОЕ H` with Сну, exit from the scattering cell апа аге 
focused by the axially symmetric non-uniform field of the focusing 
magnet to a point on the magnetic field axis, a distance Zo from 


the scattering cell. The detector which is located at 2 15 ап 


electron multiplier capable of gains up to 10 . 


А. THE DUOPLASMATRON 

The duoplasmatron is shown in Fig. 6. With hydrogen in the 
source and filament power at 35 watts, one has about .1 amp of 
electron emission from the filament. An arc voltage of 400 volts 
between the filament and the Z-electrode accelerates these elec- 
trons through the hydrogen gas, ionizing the hydrogen and forming 


a plasma. Once the plasma has been formed, an arc current of 
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Figure 6. DUOPLASMATRON 


l ampere can be sustained between the filament and the Z-electrode 
by a voltage of about 250 volts. 

To extract the plasma, the positive arc voltage is switched 
from the Z-electnodewto the anode, (see Figa 7} which is itself at 
a voltage E above ground. The energy of the hydrogen ion beam thus 
extracted will be (Ee) for singly charged ions. 

To enhance the beam intensity a simple accel-decel system is 
employed. Figure 7 shows lens l at some arbitrary negative po- 
tential V. Thus the hydrogen ions are actually extracted from 
the duoplasmatron at an energy (V+E)e. A "sock" is fitted over 
the focusing and deflector lenses and attached to lens 1 so as to 
float at this negative potential V, and thus prevent the beam from 
"seeing" ground potential in this region. Thus the beam of hydro- 
gen ions passes through the focusing and deflector lenses at an 
energy (V*E)e. Once the beam leaves the region of the sock it 
"Sees" ground potential on the walls of the vacuum chamber and 
slows to an energy Ee and passes on to the mass analyzer. 

This simple accel-decel system increased the beam intensity by 
a factor of 40. Thus we were able to obtain an analyzed beam of 
100 eV protons (with diameter .5 cm at the scattering cell 80 cm 
from the source) which was slightly greater than 1 x 1077 amps. 
Carter [6] has shown that the energy spread of the proton beam is 


less than 2.5%. 
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В. THE MASS ANALYZER 

The hydrogen ion beam is mass analyzed in the analyzer con- 
structed and calibrated by Strohsahl [7]. It employs a 45° 
bending angle (see Fig. 8) and has a resolution of n s AO wwiktb 
beam transmission of about 90%. 

Figure 9 shows the current through the mass analyzer magnet 
that will pass any particular species in the energy region from 
50 to 650 volts. Hence .45 amps allows a 100 eV beam of protons 
to travel into the scattering cell, whereas the н; апа н. аге пої 
bent enough by the magnetic field to be transmitted (see Fig. 5). 

The mass analyzer is carefully aligned with the scattering 


cell and the detector, all of which are positioned on the axis 


of the focusing magnet. (See section F below. ) 


C. THE SCATTERING CELL 

The scattering cell, shown in Fig. 10, was designed and 
tested by Bush [1]. The target gas, methane, is bled into the 
scattering cell and is accurately controlled by a Variable Leak 
valve constructed by Granville Philips. The pressure of the m 


3 


gas in the scattering chamber (approximately 10 ~ torr) is mea- 
sured by the VG ЈА ion gauge. The gas is of research grade, 
(purity 99.65%). 

The incident proton beam passes first through an adjustable 
aperature (generally set at .5 cm diameter) and then into the 
scattering chamber. Неге the H` reacts with the Е to form н, 


whose escape from the scattering cell is limited by the geometry 


of the exit slit. The exit slit is an annular opening of 3607, 
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reduced by support legs to about оон around the scattering cell, 
Separating the front from the rear of the scattering cell. The 
exit wall of the front section of the Sseatteming cell Slopes up- 
wards at 49° and*the exit wall of the reamsectin of the scat- 
tering cell slopes upward at 36". Appendix III of Ref. 5 shows 
how these two angles determine the target thickness. By varying 
the separation of the front and rear sections of the scattering 
cell, the target thickness can be increased from O to 0.8 cm. 

The 10ns collected by the front and rear sections of the 
scattering cell are measured by separate ammeters. The beam col- 
lector collects the unscattered beam that passes through the 
scattering chamber. This current, plus the current from the rear 
section of the scattering cell, is recorded as the total proton 
current incident on CHG CH), target molecules. 

Three grids on the beam collector are given potentials to 
Suppress secondary electron emission and prevents slow background 
ions from reaching the collector plate. The beam collector and 
its grids are placed inside a metal cup which shields them from 
other charged particles and also acts as a gas cap over the rear 
of the scattering chamber. 

The beam collector assembly 1S mounted on a rod which passes to 
the outside of the vacuum system, allowing the assembly to be 
pulled away from the rear section of the scattering cell and swung 
out of the path of the incident ion beam. This allows alignment 
of the entire apparatus and also allows the incident beam to fall 


directly on the detector when desired for calibration purposes. 


ДО 


р. THE FOCUSING MAGNET 

The magnetic field of the focusing magnet 1S nonuniform and 
axially symmetric and has the focusing property of bending the tra- 
jectories of charged particles so they follow a path similar to 
phat shown in Fig. lI. 

If one considers the simpler case of an axially symmetric 
uniform magnetic field, it can be shown that the trajectory of a 
charged particle is a helix (see Ref. 5). The particle will cross 
the field line on which it originated after one complete revolu- 
tion. Then, as shown in Fig. 11, particles of charge q leaving 
the source with momentum p = mv at an angle 9 will cross the 


Magnetic axis at a distance Zo given by: 


- стер COS O 
Ze qB E 


where B is the magnetic field strength. Thus a detector located 
at Z Will only sense these ions scattered out of the scattering 
cell placed at Z - O which have momentum p, charge q, and scat- 
tering angle 9. If one now holds the detector fixed at Zo and 
increases B then particles of larger p cos O values will strike 
the detector. Alternately by holding B fixed and varying 2 ме 
can detect particles scattered from the source at different values 
of p cos Ө. Since equation 1 is independent of azimuthal angle 0, 
the entire 360° can be observed at Z For most conventional scat- 
tering experiments 40 ~ "bl Hence, the solid angle is enhanced by 
a factor of approximately 300. 

In our experimental apparatus the magnetic field is NOT uni- 


form, as it is produced by a "thin" solenoid shown in Fig. 5. 
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This complicates the trajectory equations and makes numerical 


intergrations necessary. The results of numerical intergrations 
d. 3 : ‘ 
fer the H, ion are shown in Fig. 12. However, equation 1 the 


result obtained for the uniform field case is qualitatively the 
same as in our NON-uniform field configuration. 

The magnetic field of the focusing magnet was measured by 
Kelly E Figure 13 shows the axial component of the magnetic 
field on the magnetic axis. Figure l4 shows the values of the 
axial component of the magnetic field measured at various r and 
Ø values at the center of the coil. Figure 15a and b show the 
radial components of the magnetic field as measured with a Hall 
Probe(at Z, = 34 cm) across the geometric axis of the coil, both 
vertically and horizontally. This demonstrates that the magnetic 
and geometric axis of the magnetic coil are aligned to within 
0.1%. This allows use of the coil geometric axis and the coil 
center-line as reference for alignment of the entire system along 


the magnetic field axis. 


E. THE SCATTERED ION DETECTOR 

The beam detector and amplifier is a Model 306 Magnetic Elec- 
tron Multiplier whose schematic is shown in Fig. 16. The multi- 
plier is connected to a Bendix 1122 Power Supply as shown in 
Fig. 17. The output of the multiplier is fed to the input head 
of a Keithly 640 Vibrating Capaciter Electrometer. 

The Multiplier, developed by Goodrich and Wiley, provides re- 
producible current gains up to 102 when operated at the potentials 


shown on Fig. 17. When operating as the detector of scattered 
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Figure I7. Circuit Diagram for Electron Multiplier. 


ions in our experimental apparatus, the multiplier has to be oper- 
ated in a magnetic field of up to 200 gauss. The effect of this 
magnetic field is to reduce the gain of the multiplier; increasing 
field, causes the gain to decrease. 

The multiplier is operated inside an aluminum shield to pro- 
tect it from random background. The various strip potentials were 
chosen to give the highest gain and the lowest noise. The multi- 
plier gain was measured at different magnetic field values in the 
manner listed in Appendix I and fitted by a 12th order polynomial. 
See Fig. 18. The fitted polynomial was used in the actual cross 
section calculations. 

Bush, has shown that the gain does not vary with gas pressure 


up-to 2:0 X 107^ 


torr, a pressure much higher than the multiplier 
ever encounters when data is being taken. The gain is independant 
of the incident ion mass and of the incident current, provided 
the anode current is below its saturation value of 10-2 amperes. 
To prevent slow ions from reaching the detector, a grid is 
placed over the detector shield aperature as shown in Fig. 19. 


The grid also permits a crude analysis of the scattered beam energy. 


See section A, Chapter IV. 


Е. SYSTEM ALIGNMENT 

In order to detect particles scattered at a particular angle 
the detector must be placed on the focusing magnetic field axis at 
the point Z where the scattered particles focus on the axis. The 
program SOLANG which calculates the intercept distance assumes 


that the scattering took place on the axis at the center of the 
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Figure 19. Dector Aperature Geometry. 
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focusing magnetic field. Therefore accuracy in the angle of scat- 
ter determination depends on how well the scattering center and 
the detector are aligned to the magnetic field. 

Figure 15 shows that the geometric center line of the focusing 
magnetic coil is collinear with the magnetic field axis to within 
O.1%. This geometric centerline was used as the reference for 
aligning all components to the magnetic field axis. To perform 
the alignment, cross hairs were placed on either side of the mag- 
netic spool to mark its axis. Then a small laboratory laser was 
positioned so its beam lay on the magnetic axis established by the 
two crosshairs. The Br were then removed and the line 
established by the laser beam was used as the reference for 
aligning all other components to the focusing magnet axis. This 
method made possible the alignment of the scattering cell, the 
detector, and the focusing lenses to the magnetic field axis. The 
mass analyzer was constructed by Strohsahl, (see Ref. 7) so that 
the ion beam, after analysis, travels down the geometric axis. 
Hence attaching crosshairs to the front and rear of the mass 


analyzer, it was adjusted into alignment with the laser beam. 
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IV. THE EQUATIONS DESCRIBING THE MOTION OF THE 


E 
Н, ION AND THEIR SOLUTION 


А. THE TRAJECTORY EQUATION 


“Б 


The equation of motion of the H, 


ion in the nonuniform axially 
symmetric magnetic field of the focusing magnet, (equation (5)), 
is expressed in terms of the magnetic vector potential A. There- 
fore the measured magnetic field B of the focusing magnet must be 
related to this quantity before the trajectory equation can be 


solved. The expansion of A in terms of the field B is shown by 


Bush [5] to be: 


. отты à^"B (rz0,z) 
Age c cub: CARE xl (1) 
д. (көз) 0 (2) 
A (x,2) = 0 (3) 


Gagliano [9] fitted the actual magnetic field measurements along 


the r 2*0 axü with a Path "er@eér olynomial in z such that: 
I 1$ п-1 
=0 = — 
B (x A 10 Uum C2 (4) 
where I is the current through the magnetic coil. Hence we know 


В (т = 2) апа Ар ( 


Gagliano, and subsequently, Bush, has shown that the equation of 


r,z) for anysguvengE6osusung magnet curient I: 


motion of a particle of charge q and mass m with momentum p in our 


field configuration is given by: 
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dor Л 
2 дА дА 
dz ‚ 2 2 dr Ø Ø 
mv СК - Ад) - а | № 5 |* №5 ^° (5) 
| 
dz 
2 2 
where k = = (6) 
а 


To solve this second order nonlinear differential equation by 
means of the digital computer, it must be reexpressed as two first 
order differential equations. This is done by rewriting equation 


(5) as follows: 








а> (а> DNA МЫ айы ТЫП 22 (7) 
g g 
ЗА ЗА 
v (8), 5] 
Gor dz Ø ar 
__ 2» јЕ | 2 2 
k -Ag к - d 
Define: BERY (1) = г 
рввү (2) = 5% 
ЗАС 
Ag д2 
FUNC (4) = 5 5 
k s Ag 
FUNC (5) = 2 ar 2 
k e Ag 
Then equation (7) becomes 
-—— | DERY(1) | = DERY(2) (8) 


ЭЭ 


- [ppRy(2)] = DERY(2) x FuNCc(4) + [DERY(2)]°x FUNC(4) (9) 
- FUNC(5) - [peRv(2)] ^x FUNC(5) 


We now have two first order nonlinear differential equations 
which can be solved by the use of the digital computer. 

From equation (6) we see we need to know the momentum of the 
scattered particle. But p =\| 2mB_ where Е. is the energy of the 


scattered particle. Hence 


= === t 
k E! (10) 


where Е. is the energy of the scattered particle expressed in 
electron мое лао it Е. is known, equations (8) and (9) can 


be solved to yield r as a function of Z. 


+ 
B. THE ENERGY OF THE Н, 


: ‘ 


) | | + 
To determine E the energy of the scattered H 


2 particle, 


consider the following sketch of the kinematics of the reaction: 


E 
H + CH, ^ H, + СН 


Before Collision After Collision 


E : : ; 
where пы represents the mass of theH , the incident ion, and the 


m. represents the mass of the Sm the target molecule, and m 


= 


represents the mass of the Н; 


the scattered ion, and the m 
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ат 
2 А дА 
dz 4 8 Ø x 
t3 у“ Ag! Ад 7аг |“ Ag dz 2 (5) 
| + ј === 
dz 
2 2 
where k = (6) 
q 


To solve this second order nonlinear differential equation by 
means of the digital computer, it must be reexpressed as two first 
order differential equations. This is done by rewriting equation 


(5) as follows: 





slay ace] | bw et 
da (а) dz Tm № ‚ \92 [Ар 52 7 
2145) 2” А2 “= (7) 
0 g 
ЗА 9A 
В 
Agar _ @ е 
а а 
k“ -Ag 
Define: DEREI Sr 
DERY (2) = с 
дА, 
Ag Sz. 
FUNC (4) = 5 5 
к%- йн 
T 
Я дг 
Ше” с у= = 
к - Аб 
Then equation (7) becomes 
=  гввү(1)| = DERY(2) (8) 
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=- [peRY(2) | = DERY(2) x FUNC(4) + [ DERY(2) | ^x FUNC(4) (9) 
Ес 5) [DERY (2)] ^x FUNC(5) 


We now have two first ov ech ten | tats differential equations 
which can be solved by the use of the digital computer. 

From equation (6) we see we need to know the momentum of the 
scattered particle. But p =\| 2mE_ where E is the energy of the 


scattered particle. Hence 


к“ = — gp! 
В! (10) 


where Et is the energy of the scattered particle expressed in 
electron volts. Thus if Е. is known, equations (8) and (9) can 


be solved to yield Bas a function ОТ 7. 


+ 
В. ТНЕ ЕМЕКСҮ ОЕ ТНЕ H, 
| ` $ b: 
To determine ЕС, the energy of, the scattered H, particle, 


consider the following sketch of the kinematics of the reaction: 


+ 
O но Җан 


Before Collision Afer Collision 


т ь ў Е 
where м. represents the mass of the H , the incident ion, and the 


m. represents the mass of the E the target molecule, and m. 


=> 


represents the mass of the Н, 


the scattered ion, and the m 
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represents the mass of the СН, the remaining molecule. Ignoring 
the binding energies, we can write the following: 


From the conservation of energy we have; 


п о о m v (11) 


From conservation of momentum we have: 


пту = ту соз Ө + шщ v cos f (12) 
о 5 5 ЕЕ 


о а Е ; 
m v. sin 0 - m,v, sin Ø (Л) 


We wish to solve these three equations to yield the relation: 


Squaring (12) and (13) and adding, we eliminate Ø and get: 


22 DM p 2 
mv - 2m vmv cos ОӨғптпу = му 
O о O O S S S S r r 
2 2 2 2 
But from (11) m v =ппу -тту . 
т гоо rs s 


Thus we have on substituting this into the above equation: 


2s -mm + пт й2 = 0 
$ $ тео © 


DUE 
mv - 2m vmv сос Ө +m 
оо 5 ír ss 


О О 5 


or 2m E  - Б Em E сов @ + 2m E_ - 2m E_ + 2m_E_ = 0 
— оо оо $ $ 8 5 r o r s 
Thus we get on rearranging: 


E (m, * т.) 4; E (m, - т.) = 2 ee cos 9 (14) 


Now let: а= т * m, and b = mo - m, (15) 


Squaring (14) we get: 
2 2 WD 
A EA + |2ab Е - & юш ш Е cos ој E +b E = О 
S о ое о 5 о 


, 


p 


Solving for E we obtain: 


2 Ше 2:22 02 
|l mm E соѕ Ө - 2аЪЕ + \{(2аЪЕ - 4m_m_E_ cos 9) - la b E 
= O S о o — О o so О 
2 
7 2а 


As the + root only yields physical results, in our case we have, 
on substituting equation (15) for a and b: 


2 
2mm_cos © + (m_-m_}(m_+ m_} 
Е о 5 т шо - r 
E = Е |1-------------------- + 
Б О 2 
(m. * m) 
S E 


A m^n^ cos tome Ne cuc) TURIN со 0 
"o s s "r/ V 99m os 
2 
Ez 
(m + m) 
Now: m, = ag 
2 
m =m 
О HË 
= = 2 
m. m m 
2 
Thus the above equation simplifies to: 
2 
n^ * m m A 2m + hm cos Ө 
= --------------------- 16 
в А (10) 


(m, + 2m ) 


\ 2 2 
2 |2 м соѕ Ө me + mm - Be + 2m cos Q 
| О r о О О 





= 
2 
+ 
(m, 2m.) 
Thus selecting 9 and E we can calculate Ес, the energy of the 
T 5 
scattered H,, from equation (16). 
Е 
Figure 20, a plot of the ratio = versus the scattering angle 
О 


9, shows that in the oh range from о Eile Е. changes by less 


than 3%. 
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.948 


Figure 20. Ratio of Scattered H5 
Energy to Incident 
Ht? Energy versus 
Scattering Angle. 
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С. THE METHOD OF SOLUTION 
Assume protons of energy Во and a focusing magnetic current ІЗ 
to have been selected. Then equation (4) yields the magnetic 
s . Е + 
field В (т = 0,2). Equation (16) yields the energy of the H, as 
2 


a function of the scattering angle, O. From equation (10) k is 


4- 
then known for the H, and we are ready to solve equations (8) and 


2 
(9). The solution is a trajectory similar to the ones shown in 
Fig. 7. The distance Zo where r is zero is recorded, This is 


where the detector must be located to detect н, particles scat- 
tered out of the scattering cell at the angle 9 when the focusing 
magnet current is set at I, and the energy of the incident beam 
of protons is EG: 

If the detector distance Zo 1s changed and 16 апа E. are held 
constant then H, particles scattered out of the scattering cell 
at different scattering angles 9 аге detected. 


Alternately, one can hold Eo and Zo constant and vary IS. This 


also allows particles of different scattering angles to be detected. 


D. THE COMPUTER SOLUTION 

The computer solution of equations (8) and (9) was originally 
worked out by Gagliano, but has since been rewritten. Program 
SOLANG, explained in Appendix II, employs the Hamming Modified 
Predictor-Corrector method to integrate equations (8) and (9). 
This is the DHPGC subroutine of the IBM System 360-67 Scientific 


Subroutine Package. 
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The trajectory equations are solved for a range of scattering 
angles for a particular Ey and Io: The results are stored and 
later used to compute the solid angle subtended by the detector as 


it intercepts each of these trajectories. 


Е. THE SOLID ANGLE 


The solid angle subtended by the detector is given by the 


relation: 
а0 (0,0) = sin O dO dağ 


But since our system is axially symmetric this becomes 


аб (о) = sin 0 do fag 


- 2m (225) sin 0 do 


1.612 " sin 0 40. 


This is because only 290° of the available 360° of the axial angle 
Ü are accessible to the detector due to the construction of the 
scattering cell. (See Chapter II, Section 3.) 

From the trajectory plots on Fig. 15 it can be seen that all 
particles scattered into 9 will cross the axis at z within AZ. 
This AZ determines the angular resolution AQ which the detector 
Sees and so allows determination of the solid angle. But AZ is 
in turn determined by the detector geometry and the particular 
trajectory under observation. Тһе relation between ЛӨ апа AZ can 


be found from the functional relationship: 


Zo = f(E, B, о,0 2 a, м 


where E = energy of incident ion. 
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B = magnetic field 

“= target width from system axis along a radial line 

N = target thickness from system center along the Z axis 

C = constant for each reaction including mass ratios, 
inelastic energy losses, etc. 


Differentiating this gives: 


O 


92 92. д2 02. 92, az. 
4. шы "ы ӨЗ аа Алани". tye АС 


But AB = O and AC = O for any specific reaction and measurement. 
Thus, solving for AQ we get: 
à ға) | д Žo) | in 
MM Ru пы Е (= ИР = AT) 
o \ОЕ д 
AQ = 0,0,1 Ё Е,0,1 91 Е,0, 
92 


Е 


No analytical expression exists for the various differentials. 


ES) 


However, computer trajectory results obtained by varying HEMOS 
and T] show the effect of the differentials is to change AO by no 


more than by 10% (see Ref. 10). Thus the angular spread can be 


AQ ~ (85) AZ. [51477 
E 


O 


expressed as: 


Hence the solid angle becomes 
: o9 
dO (o) z T OPI sun Qf - e fs) 
92 6 о 


Мои, ѕіпсе AZ depends on the angle РБ (see Fig. 21 and 
Ref. 5 for further details) and 9. in turn depends on the scat- 


tering angle 9 we see that d{!must be evaluated at each angle of 
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Figure 21. Angular and Energy Acceptance of Dector 
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scatter. This is implemented in the latter subroutine of the 
computer program SOLANG. 
Computer results indicate that the angular resolution given 


by equation (17) is approximately Ли 


F.— THE CROSST SECTION 
Consider a beam of monoenergetic ions I; incident ona target 
of Ny molecules. A detector located at an angle 90 subtends solid 


angle dQ (0). 





The number of particles scattered into the detector, ve? 1$ 


usually measured as a current I. where: 


І. = м 40 (0) 


The cross Section is then defined as: 


I 


5 
(9) * Y-N. dü(6) (12 
T t 
But I. is magnified G times by the multiplier gain. Hence Ip: 


the actual detector current measured is: 
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Hence (19) becomes: 


I; G N, а0(0) 


But io the number of target particles, is given by: 


Where р = pressure of target gas in scattering cell in torr 
16 
= 3.586 х 10 Бап one torr 


t = target length in cm 


Hence we have: 


1% 2 


—————— M cm (20) 
Іі. СР 0% аП(0) 


do (0) = 


Because of the construction of the scattering cell (see 
Section 3 of Chapter II) the target thickness depends on the 


Scattering angle O. Bush, has shown this relation to be: 


О 
ме 2.54 <іп 136 -0 if o < 36° 


Sin 9 
ї = $ if 369 « 0 « 49° (21) 


, о 
=e 2.0 sin (9-49 if 99 С 


Sin 0 
where S is the minimum separation of the front of the scattering 
cell from the rear of the scattering cell, (which is adjustable 
from outside the vacuum system). 
Computer Program CSVSZ, listed in Appendix I, calculates the 
differential cross section in the Lab. and Center of Mass coordin- 
ate systems using equations (20) and (21) and the solid angle 


determined by the program SOLANG. 
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У. EXPERIMENTAL RESULTS 


А. ELIMINATION OF BACKGROUND 

When н” collides with CH, three general types of collision may 
occur which can produce large detector signals. 
Event (a) The rearrangement collision where the н" captures a 


+ 


proton to form H, 


Event (b) A scattering collision of the H with a hydrogen atom. 
Event (c) A scattering collision of the H' with a carbon atom. 
Events (b) and (c) are expected to have cross sections that are 
many orders of magnitude larger than those of event (a). Hence 
appropriate steps were taken to insure events (b) and (c) were 
not measured by our detector. То eliminate events (b) and (c) 
we look first at the energy and momentum of the scattered parti- 
cles in all three events. 

Consider first event (a). The capture process may be approxi- 
mated by the reaction shown below: 


+ т 
ы) а Сна јен 


= p 
V 
S 
mo ж. 
M \ зи СР 45° 
О а... ЕЗ 
mo Do "x 
E: 
И 5 + А : 
The energy of the incident H particle is: 
2 
2 Po 
- 1 А 
оо = 2m) (1) 
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Now, from conservation of momentum we see the momentum of the 


+ 
scattered H, particle, represented by p.» related to M by: 


P 
о Б 
Ро =Р. cos 45° = = (2) 


where we assume the mass of the CH represented by m. , to be 


3? 


infinitely large. Hence in the capture process, the scattered 


particle has momentum ES 2 Ро. The energy of the scattered 








particle is then: 2 2 
E = Ps Ро 
5 2m ~ m 
S 5 
+ : 
Büt for the H. particle; m = 2 mo 
2 
Hence: Бе 
в. = 2m = Во . (3) 


i ; d 
Consider now event (b), the scattering of the H from a 


hydrogen atom: 


V 
S 
m о 
Ма 5 РА Ө ~ &5 
tg "Осы mer 
О 
m. S 
t 
Pt 
Since mn, = м. апа р, = р, ме ѕее that: 
о 
Ро = 2p, cos &5 - (2 P. 
Po 
Hence: р. = = (4) 
2 2 
Тһеп: и о КЕ #5) 
| 5 2m. Ат PON 


3 s : 
Thus the elastically scattered H from H has a momentum that is 


1 
Е 


approximately —^ and an energy that is approximately Е 


2 
. а + 
Finally consider event (c), the scattering of the H froma 


carbon atom: 


+ + 
H +C>H +C 





Assuming the carbon to be infinitely massive we see that: 








Е. = E (6) 
2 2 
Непсе: Рз " Po 
2m 2 
S О 
But: m =m 
5 О 
Thus: BaS Ro (7) 


Hence the elastically scattered H from C has momentum Po and 
energy Ес: 

Summarizing the results ОҒ this very approximate analysis in 
Table 1 we get for an incident beam of protons of energy Е апа 
momentum Ро the information shown. Note that the scattered H' 
produced in events (b) and (c) has momentum less than the н> 


а 


produced in event (a). 
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Possible Events Scattered Momentum of Energy of 
Расе Scattered Particle} Scattered 
Particle 


Event (a) Capture 
E ишен. * CH н? 12 Е 
а & ^ E Fe 


Event (b) Scatter 





+ + + " 
Рон њи H “в. 
Event (с) Scatter 

+ E + 
Н “Сен «Сс H р E 
= = о О 


TABLE 1. The Three Possible Scattering Processess. 


Plots of the scattered ion trajectories resulting from each 
of the three events tabulated are shown in Fig. 22, for the parti- 
cular case of E. = 100 eV and I = 6.0 amps. Note that the two 
elastically scattered H“ ions cross the axis at a Z less than that 
of the H, ion. To prevent these unwanted H' particles from 
reaching the detector, a baffle, 5 cm in radius was positioned 
over the beam collector as shown in Fig. 22. This allowed only 
the H, produced by the rearrangement collision to reach the 
detector. 

Figure 23 contains plots of the detector current versus the 
detector grid voltage for the case of E - 100 eV and I - 6.0 amps 
and with the detector located at Z = 46 ст (this correspemds to a 
scattering angle 9 of 46.927). Figure 23 -(i) was recorded before 
the baffle was placed in the path of the E fens and Fig. 23 - (11) 


Was recorded after the baffle was positioned, as shown in Fig. 22. 
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Consider først, Figo29 1| ОТТА асос ОЛЕ 
eliminates the slow ions. As the detector grid voltage is varied 
from 25 to 50 volts, we continue to collect Hhewmons Sea titer cama on 
the scattering cell due to processess (a) and (b) and (c). Between 
60 and 90 volts we collect only these particles with energies 
greater than 90 volts; ie, the ions which result from processes 
(a) and (са ошу 

Figure 23-(ii) now shows that the baffle does indeed eliminate 
particles resulting from processes (b) and (c). Thus with the 
grid voltage set anywhere between 30 and 90 volts, and the baffel 

+ 


in position, we only detect the н, that results from the re- 


arrangement collision. 


B. DEMONSTRATION OF VALIDITY OF DETECTOR SIGNAL 
From equation (20) of Chapter IV we have: 


in 2 


a ЧЕ (8) 
I; P... 616 99(0) 


ас(8у- 


If the angle of the scatter O is held constant, then this re- 


duces to: 


D 
с = ------ 
а KA n 
i sc 
where K - constant. 
Consider now the case where І; is held constant. Then (9) 
reduces to: 
= Б 10 
1p n Sc ( ) 
where C, = constant. Experimental data plotted in Fig. 24 obeys 


1. 


this linear relationship. The fact that I, does not extrapolate to 


zero 1S due to the background. 


70 


|9) Би!2э4405$ 
(шо) 7 2042909 шоәа 


46 %6 16 8, Se 25 66 96 ce ос ¿Z z 12 ШЕТ 
Ns 


>< 2143 08 9. 
Г 


4- E NN 6 
LULA LL LL LL LLL AL LALLA LLL LLL ALLTEL I LL LL LU ULL LU LLLI IL ULL IL ULL LU LL UL I UL LLLI LLL LLL LU LLL UU di 777777777777 


” A 


лэашоцо шппэол 4O SIOM 


с О 


9 
WIA 
€ 









we 
(9) мәлз| .Н (а) eng] oH 


(о) мәлз| ён 


а 


(шә) ð 
лә 001- шрад цоуола уцәр:ди! 


ду, јо Аблаџа ду рир 'дашо Q 9 - | чәчм 
oGt - 10 S9|9!JJDd ,H p9349jjD9S 5||021150|Ә 


au, puD ?H 8uj jo sjoid K10j2efDJ] *ez aJnDiJ 


Detector current (Amps) 


io * 


Под 


Тар 


%0 


Figure 23. Detector current vs. Detector grid voltage 
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Consider next the case where бы is Weld constant. Then 
equation (9) reduces to: 


pe -wc m (11) 


where 1. = constant. Experimental data plotted in Fig. 25 follows 
this linear relationship. Thus our detector signal has the proper 
dependance on the various parameters and behaves like a scattered 
particle current. 

For the differential scattering measurements made with 100 еу 
incident protons, the proton current at the beam collector was 
approximately 107” amps and the current to the rear of the scat- 
tering cell was about о" amps. The methane gas pressure was 


4 3 


ЧЧ еа from 1 x 10 ' torr to 3 x IO " senmr. The ааа signal 
ONG -12 
under these operating conditions varied from 1 x 10 amps to 
-11 
вк 10 amps. This represents a flux of a few hundred particles 
per second. With no gas in the scattering chamber the detector 
ES 

Current fell to about 5.0 х 10 amps. 

The UGIA Ion Gauge used to measure the methane pressure in the 


scattering cell was calibrated using a Capacitance Manometer. The 


relationship between the pressure measured by the ion gauge 


Ge ) and the actual pressure which was measured by the mano- 
measured 
meter (P ) is linear in the region of interest as is shown 
actual 
below. 
m +" 2 slope = 0.625 
Э 


Па холо" 
t OD 


РО 2.0 5149, 


P x 19 Тоту 
measured 


ps 


атр 


12 
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Figure 25. Detector Current versus Incident Current. 
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С. ss CH, SCATTERING DATA 

The cross section for the Тоғтастон оя H, was measured at 
different angles of scatter by varying the detector distance Z 
while holding the magnet current I fixed. As pointed out in 
Section B of Chapter II the cross section is expected to be a pro- 
nounced peak at 46.9°. Hence the detector was swept through the 
range of scattering angles from nos to ho. This was repeated for 
various target gas pressures for each value of the energy of the 
incident proton beam. The energies investigated were 70, 85, 100, 
150 and 200 electron volts. The scattering data recorded at each 
of these energies is shown in Figs. 26, 28, 30, 32 and 33 
respectively. A prominent peak exists in the 70, 85, and 100 eV 
data. The 150 and the 200 eV data shows no apparent peak in the 
cross section. 

The data was analyzed in the following manner. From the cross 
section measurements made at each energy an average value of o(90) 
was computed and sketched in as a solid line. The assumed back- 
ground was then drawn as a dotted line. Subtracting out the back- 
ground from the 70, 85, and 100 eV data, the points connected by 
the solid curve in Figs. 27, 29, and 31 were obtained. To find 
out what the actual peak in the scattering data looked like before 
it was detected by our detector of approximately 1° resolution, 
the effect of the angular width of the detector was determined and 
the actual curve unfolded. The curves dotted on Figs. 27, 29, and 


31 represent the unfolded original peak in the scattering data. 
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Figure 33. Cross section for 
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47 


48 


included) 


49 


These peaks are then the actual cross sections for the rearrange- 
ment reaction: 
НЫ БУНЕ Ны ен 
= Сое 
"s 3 
as the energy of the H varies from 70 to 100 ev. All the peaks 
exhibit the same characteristics: a sharp leading edge at the lower 
angles, a flat top and a long decreasing tail extending out to 


about Б. The peaks dependence on energy is exhibited in the 


following table: 
















@ Q 
> > 
H H 
5 ч з 
~ ic: ы и о 
Q © m © Фф 
>) Q o © © Sg) 52 % 
Ер а м ч E coo ч чә 
р б 5 = ruv, e» (0) O O 
ч O Фо d д Hori 
CE Е О ~ б Qo а о, ге 
С ао с Ces QE С Е 
ej H E Q M = О “о 
(x) > О < n. < < OU O Q 
| | 
= i 
EE c E 10.2 |20.6 45.879 
-2 2 
85 eV| 4.26x10 c 6.56 11. 50 45.6” 
-22 2 О 
100 eV| 3.0x10 cm LIS 746 ДО 
= 2 ee 
2150 еу|=1.0 xl10 em Not Not Not 


Observed|Observed|Observed|Observed|Observed 


Table II. Summary of Experimental Data 


The peak does exhibit the expected physical properties. As the 
energy increases the magnitude of the peak decreases, the width 


о 
narrows and the angle at which the peak occurs moves towards 46.9 . 
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The energy dependence of the experimental cross section is compared 
to the theoretical cross section in Fig. 34. Curve (a) represents 
s ie. the magnitude of the flat top. Curve (b) represents 


the area under the flat top. Curve (c) represents the total area 


under the curve. 
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Figure 34. Energy dependence of 


cross section 


Bates, Cook, Smith 


present results 
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У. CONCEDSION 


The close agreement between the experimental data and the 
predictions of the ion molecule rearrangement theory proposed by 
Bates, Cook and Smith is evident in the following observations: 
(1) The sharp peak in the cross section predicted by Bates et.al., 
is deserved in the 70, 95, апа 109 eV data. Іре Тас И ECIS 
peak is not observed on the 150 and 200 eV data is not sur- 
prising for extrapolation of the experimental data to 150 eV shows 
that the magnitude of the cross section should be 1.0 x ПО су, 
The 150 eV experimental data on Fig. 32 shows that a peak of this 
magnitude or smaller cannot be observed. 

(2) The theory predicted the peak to occur at 46.9°. The peaks 
observed were located between и апа 46.257, the latter cor- 
responding to the higher energy. As the accuracy of the approxi- 
mations made in the theory increases with higher energy, one 
expects the location of the peak in the experimental data to ap- 
proach 46.9° as the energy is increased, This trend is observed 
experimentally. 

(3) As is evident from Fig. 34, the experimental cross section has 
an energy dependence similar to that predicted by the theory. 

(4) The cross section predicted by Bates et.al., is an upper limit. 
Hence it is expected that the experimental cross section might be 
smaller than the theoretical cross section. The experimental cross 
Section shown in Fig. 34 is in fact about a factor of 30 smaller 


than the theoretical upper limit. 
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(5) As some of the approximations made in the theory may no longer 
be valid at low energies one expects the peak to broaden as the 
energy decreases. This trend is indeed observed in the 100, 85, 
апа 70 ем data shown in Figs. 31, 29, and 27 respectively. 
It is felt that the facts mentioned above constitute a satis- 
factory verification of the ion molecule rearrangement theory of 
+ 


Bates, Cook and Smith as applied to the formation of Н, in the 


reaction 


o 


HO CH ^ H. + CH 
a 47 82 7 сң; 
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АРРЕМОТХ Т 


MULTIPLIER GAIN MEASUREMENT 
The following sketch shows how in principle, the gain measure- 


ment can be made. 


t 
t 
1.” Dynode 
JL t E Str ip 
Field 
Strip ~ Аподе 
E 


By recording the incident current Ii and the detector output Ip 


the gain can be computed from: 


_ - 
Iw. 
1 
To measure I. the Dynode Strip was used as a Faraday Cup. Two 


problems then arose: 
: 8 

(1) The Dynode Strip, because of its large resistance (10 00) 
retains a large residual negative charge. Thus in using the Dynode 

-10 
Strap as a Faraday Cup, І. had to be greater than 10 amps so 
as to insure that this residual negative charge did not adversely 
affect the measurement. 

— “6 
(2) The Multiplier saturates when In 2 5x 10 amps. Thus we must 
-6 

have І; С. 

То overcome these two difficulties the gain measurement was 
performed in the following sequence of steps. 


-10 
Step 1. A beam of about 10 1 amps was focused on the Beam 


Collector. 
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Step 2. 


Step 3. 


Step 4. 


Step 5. 


Step yo. 


The Beam Collector was swung away from the back of the 
Scattering Cell. ТАЕ оса the ion beam to fall 
directly on the multiplier. With F.S.I. and D.S.O. 
grounded, (see Fig. 17) and a Keithley 610 Micro-Micro- 
ammeter connected to the D.S.I. this incident current 

was measured to be m (= ъб 25 amps). 

The Multiplier was now restored to the configuration shown 
in Fig. 17. The Dynode Voltage was reduced to approxi- 
mately 1,000 volts (this insured that the anode current to 


. -6 
be measured in Step 4 was less than 10 ^ amps the satu- 


ration current of the Multiplier). 
1 


The Multiplier was turned on and a detector current Ip 


57 


(~ 10 amps) was measured at the anode. Thus, with the 


Multiplier settings as selected in Step 3 the gain was 


ї . 

The incident beam І. was now reduced by adjusting the 
focusing electrodes (see Fig. 5) a significant amount to 
1 . 

cause 1, to decrease at least two orders of magnitude. 


2 1" -9 
This new detector current Ip (= 10 amps) was then re- 
I 
. е р 
corded. Thus the incident beam must have been І. = G^ 5 


Тһе Dynode Voltage was turned up to 1,500 volts. Тһе 
Multiplier now had the configuration desired for the 
experimental measurements to be made later on. The de- 
tector signal Ip now recorded allowed the gain of the 


Multiplier in this configuration to be computed from: 


axes 
I. 
L 
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The above procedure was repeated at various values of detector 
distance Z and focusing magnet current I selected so as to give a 
measure Of G versus the axial magnetic field B. The results of 


such measurements is shown in Fig. 18. 
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APPENDIX ITI 


The two computer programs SOLANG and CROSEC are listed in the 
following chapter. 

SOLANG, is the program which integrates the equation of the 
motion of the scattered H, ion (equations (8 and(9)of Chapter IV) 
for predetermined values of the Focusing Magnet Current I, and 
the Incident Proton Energy E over a range of scattering angles 0. 
The integration is preformed using the D.H.P.G.C. Subroutine of 
the IBM System 360-67 Scientific Subroutine Package in the MAIN 
portion of the program, and in the two Subroutines FCT and OUTP. 
The resulting distances Z where the trajectory crosses the 
magnetic axis and the corresponding angles pe (see Fig. 21) are 
stored and later tabulated. 

The program now continues on to the Subroutine DOMEGA and the 
succeeding subroutines, to compute the solid angle d2(@) for each 
of the н? trajectories. The solid angle is computed using 


2 


equation (18) of Chapter IV, which is: 


dQ(@) = 1.6127 sin • (82) 42, 
oz 
OJ EST 


AZ. has been determined by Bush, (Ref. 5) and is given by the 
following relations: (see Fig. (19) for explanation of symbols 
used) 


I fore NES о> 26.6° 


|| 


AZ ( ftad cotan om - 2b 


92 


ТЕ fir 2696? 2 Өк 2 а 
AZ = (Һғаһт-сосбатя6” - 6 
О О 
ШЕТ бот 47.9 7%- 09 
AZ = 2h cotamnmwo 
О О 


: од : — 
The relation ($27 \ 1S determined by fitting the results of Zo 
O/E,I 


versus 9 that were obtained in the trajectory integrations by a 


OZ 


polynomial 7 с ЈЕ Е(0). We then evaluate O° at the various G 
values. In this manner the solid angle dQ is computed for various 
values of 0. 

In order that the results may be tabulated for various Z values 
it is necessary to fit the results of O versus Zo that were ођ- 
tained in the trajectory integrations by a polynomial 9 = f(Z). 
Hence, on choosing a Z, the corresponding 9 can be computed and 
in the manner described above, the solid angle can be determined. 

CROSEC is the program which computes the cross section from 
the experimental measurements of the Detector Current I5 and the 
Bmedeent Prot&en Current I; which were recorded at various detector 


distances Z for a particular proton energy E, methane pressure бе. 


1 C 


and a particular Focusing Magnet Current B 

Because the data is recorded at various detector distances Z, 
the coefficients of the polynomial O - f(Z) that were determined in 
the program SOLANG are read into CROSEC on data cards and so allow 
the O corresponding to each Z to be determined. Rather than read 
into this program the solid angle at each point Z, the coefficients 


of the polynomial А = £(9) that were determined in SOLANG are read 


in on data cards. Hence using: d(9@) = 1.61217 sin 222) AZ. and 
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the relations for AZ referred to in the discussion of the program 
SOLANG, we compute the solid angle at the O corresponding to the Z 


where the data was recorded. Hence using equation (20) of Chapter 


Ту lp 
4269) = aie 
Т.Р сс et G dQ(@) 


we compute the cross sections. 


Using the transformation from the Laboratory to the Center of 


+ i 
Mass coordinates for the reaction lieet m H, + Cale, we can de- 


j о (9 . 
termine Ө РЕ" and hence o( 1% of M The results are both 


tabulated and plotted. 
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